ON EXPONENTIAL CONVERGENCE OF GEGENBAUER 
INTERPOLATION AND SPECTRAL DIFFERENTIATION 
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Abstract. This paper is devoted to a rigorous analysis of exponential convergence of 
polynomial interpolation and spectral differentiation based on the Gegenbauer-Gauss and 
Gegenbauer-Gauss-Lobatto points, when the underlying function is analytic on and within 
an ellipse. Sharp error estimates in the maximum norm are derived. 



1. Introduction 

Perhaps the most significant advantage of the spectral method is its high-order of accuracy. 
The typical convergence rate of the spectral method is 0(n~ m ) for every m, provided that 
the underlying function is sufficiently smooth [351 111 13 HH1 IE]- If the function is suitably 
analytic, the expected rate is 0(q n ) with < q < 1. This is the so-called exponential conver- 
gence, which is well accepted among the community. There has been much investigation on 
exponential decay of spectral expansions of analytic functions. For instance, the justification 
for Fourier and/or Chebyshev series can be found in [551 113 [Ml HE [551 SI]- I n the seminal 
work of Gottlieb and Shu [2H [55] , on the resolution of the Gibbs phenomenon (the interested 
readers are referred to Gustafsson [2S] for a review of this significant contribution) , the expo- 
nential convergence, in the maximum norm (termed as the so-called regularization error), of 
Gegenbauer polynomial expansions was derived, when the index (denoted by A below) grows 
linearly with the degree n. Boyd [7] provided an insightful study of the "diagonal limit" (i.e., 
A = fin for some constant /3 > 0) convergence of the Gegenbauer reconstruction algorithm in 
|24j . We remark that under the assumption of analyticity in [24] . the exponential accuracy 
of Gegenbauer expansions is actually valid for fixed A (see Appendix [A] for the justification). 

It is known that the heart of a collocation/pseudospectral method is the spectral differ- 
entiation process. That is, given a set of Gauss-type points {Xj}J =0 , e.g., on [—1,1], the 
derivative values {u'(xj)} can be approximated by an exact differentiation of the polynomial 
interpolant {(I n u)' (xj)}. Such a direct differentiation technique is also called a differencing 
method in the literature (see, e.g., (39J [El [35] ) . For the first time, Tadmor [39] showed the 
exponential accuracy of differencing analytic functions on Chebyshev-Gauss points, where 
the main argument was based on analyzing the (continuous) Chebyshev coefficients and the 
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aliasing error, and where the intimate relation between Fourier and Chebyshev basis functions 
played an essential role in the analysis. Reddy and Weideman [35j took a different approach 
and improved the estimate in ■''>!) for Chebyshev differencing of functions analytic on and 
within an ellipse with foci ±1. As pointed out in [35], although the exponential convergence 
of spectral differentiation of analytic functions is appreciated and mentioned in the literature 
(see, e.g., [TH1 HOI H]), the rigorous proofs (merely for the Fourier and Chebyshev methods) 
can only be found in [3H1 [35] ■ Indeed, to the best of our knowledge, the theoretical justi- 
fication even for the Legendre method is lacking. It is worthwhile to point out that under 
the regularity condition (M): \\u^ < cM k , the super-geometric convergence of Legendre 
spectral differentiation was proved by Zhang |46j : 



oW J(u-I n u)\ Xj )\<c(^L_)-\ (1.1) 



where {a^}™_ are the Legendre-Gauss points. Similar estimate was nontrivially extended to 
the Chebyshev collocation method in [57]. The condition (M) covers a large class of functions, 
but it is even more restrictive than analyticity. On the other hand, the regularity index k 



could be infinite, while the dependence of (1.1) on k is not clear 



The main concern of this paper is to show exponential convergence, in the maximum 
norm, of Gegenbauer interpolation and spectral differentiation on Gegenbauer-Gauss and 
Gegenbauer-Gauss-Lobatto points, provided that the underlying function is analytic on and 
within a sizable ellipse. The essential argument is based on the classical Hermite's contour 



integral (see (2.19) below), and a delicate estimate of the asymptotics of the Gegenbauer poly- 



nomial on the ellipse of interest. It is important to remark that the Chebyshev polynomial 



on the ellipse takes a very simple explicit form (see, e.g., [ID] or (3.5) below), but the Gegen- 
bauer polynomial has a complicated expression. Accordingly, compared with the Chebyshev 
case in |35j . the analysis in this paper is much more involved. The Chebyshev and Legendre 
methods are commonly used in spectral approximations, but we have also witnessed renewed 
applications of the Gegenbauer (or more general Jacobi) polynomial based methods in, e.g., 
defeating Gibbs phenomenon (see, e.g., [2H0H]), hp-e\emeiits (see, e.g., [TJJ [2 EH US] ) , and 
numerical solutions of differential equations (see, e.g., [3] fE [ 127 ] fH ] fT3 ] [Ti l IT5 ] fTH ], EH ] 137] ) 
and integral equations (see, e.g., [§])■ The results in this paper might be useful for a better 
understanding of the methods and have implications in other applications. 

The rest of the paper is organized as follows. As some preliminaries, we briefly review basic 
properties of Gegenbauer polynomials, Gamma functions and analytic functions in Section 2. 
We study the asymptotics of the Gegenbauer polynomials in Section 3, and present the main 
results on exponential convergence of interpolation and spectral differentiation, together with 
some numerical results and extensions in the last section. 



2. Preliminaries 



In this section, we collect some relevant properties of Gegenbauer polynomials and assorted 
facts to be used throughout the paper. 
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2.1. Gegenbauer polynomials. The analysis heavily relies on the normalization of [3B], so 
we dehne the Gegenbauer polynomials^ by the three-term recurrence: 

nC x (x) = 2(n + \-l)xC x _ 1 (x)-(n + 2\-2)C x _ 2 (x), n>2, 

C${x) = l, C x (x) — 2\x. A>-l/2, jce[-l,l]. 

Notice that if A = 0, C„(x) vanishes identically for n > 1. This corresponds to the Chebyshev 
polynomial, and there holds 

2 2 

lim A _1 C^(x) = -TJx) = - cos(n arccos(x)), n > 1. (2.2) 
X-yO n n 

Hereafter, if not specified explicitly, we assume A 7^ 0, and refer to [35] for the analysis of the 

Chebyshev case. Notice that for A = 1/2, C x (x) = L n (x), i.e., the usual Legendre polynomial 

of degree n. 

The Gegenbauer polynomials are orthogonal with respect to the weight function (1 — 
x 2 ) x ~ 1 / 2 , namely, 

/•1 

C x (x)C^{x){l - x 2 ) x ~ l ' 2 dx = h x 8 mn , (2.3) 

where S mn is the Kronecker symbol, and 

fe A_ ^r(n + 2A) 
" r 2 (A) n!(n + A)' 1 ' 

Moreover, we have 

c x (-x) = {-i) n c x (x), gn(i) = r ^r| 2 Ar (2 - 5) 

and 

^C x (x) = 2\C x ±l(x). (2.6) 
By Formula (4.7.1) and Theorems 7.32.1 and 7.33.1 of Szego [38], we have 

\C x (x)\ < C x (l) if A > 0; \C x {x)\ < Dxn*- 1 if - ~ < A < and n > 1, (2.7) 

where D\ is a positive constant independent of n. A tight bound can be found in |34j (also 
see |32j): 

max{(l-x 2 ) x (C x (x)) 2 } < 2e (2 + V2*) fe A A > , „ > 0. (2.8) 

2.2. Gamma and incomplete Gamma functions. The following properties of the Gamma 
and incomplete Gamma functions (cf. [45j ) are found useful. The Gamma function satisfies 

T{x)T{x + 1/2) = 2 1 ' 2x y/^T(2x) 1 Vx > 0, (2.9) 

and 

r(x)T(-x) = ^ — r, Vx>0. (2.10) 

x sm(7ra;) 

We would like to quote the Stirling's formula (see, e.g., [215]): 

V2nx x+1 / 2 e~ w < T(x + 1) < \f%KX x+l l 2 e- x e^ , Vx > 1. (2.11) 
We also need to use the incomplete Gamma function defined by 

/>oo 

T(a,x)= / 1?- 1 e- t dt, a > 0, x > 0, (2.12) 



^Historically, they are sometimes called "ultraspherical polynomials" (see, e.g., the footnote on Page 80 of 
[38] and Page 302 of [l]). 
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which satisfies (see P. 899 of [3D]) 

n 

r(n+l,a;) = n\e~ x ^2 



k=0 



X 



n = 0,1, 



(2.13) 



2.3. Basics of analytic functions. Suppose that u(x) is analytic on [— 1, 1]. Based on the 
notion of analytic continuation, there always exists a simple connected region R in the complex 
plane containing [—1,1] into which f(x) can be continued analytically. The analyticity may 
be characterized by the growth of the derivatives of u. More precisely, let C be a simple 
positively oriented closed contour surrounding [—1, 1] and lying in R. Then we have (see, e.g., 

my- ' 



\u^{x)\ < max ze c\u(z)\L(C) 



Vx e [-1,1], 



(2.14) 



where L(C) is the length of C, and 8 is the distance from C to [—1,1] (can be viewed as the 
distance from [—1, 1] to the nearest singularity of u in the complex plane). Mathematically, 
an appropriate contour to characterize the analyticity is the so-called Bernstein ellipse: 



E p := [z 



z = —{w + w 1 ) with w — pe l9 



e 



[0,2tt]}, 



P>1, 



(2.15) 



where C is the set of all complex numbers, and i = \f— T is the complex unit. The ellipse E p 
has the foci at ±1 and the major and minor semi-axes are, respectively 



1 { P + P - 1 ), b^lip-p- 1 ), 



2 yr r n 2 
so the sum of two axes is p. As illustrated in Figure 



(2.16) 



2.1 



p is the radius of the circle w = pe" 



that is mapped to the ellipse £ p under the conformal mapping: z = \(w 



|w|=p 





Bernstein 
ellipse 








— ^ > 







-1 







1 



-1 



1 



FIGURE 2.1. Circle (left): \w\ — p — 1.5, and Bernstein ellipse (right): £ p with 
foci ±1 linked by the conformal mapping: z — =(w + w _1 ). 

According to [3S] , the perimeter of £ p satisfies 

L{£ P ) < ttVp 2 + P- 2 , (2-17) 

which overestimates the perimeter by less than 12 percent. The distance from £ p to the 
interval [—1, 1] is 



5 P = nip' 



(2.18) 
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We are concerned with the interpolation and spectral differentiation of analytic functions 
on the Gegenbauer-Gauss-type points. Let {xj :— Xj(X, n)}™ =0 be the Gegenbauer-Gauss 
points (i.e., the zeros of C* +1 (x)) or the Gegenbauer-Gauss-Lobatto points (i.e., the zeros of 
(l-x 2 )C*±{(x)). The associated La grange interpolation polynomial of u is given by I n u e P n 
(the set of all polynomials of degree < n) such that (I n u){xj) = u{xj) for < j < n. Our 
starting point is the Hermite's contour integral (see, e.g., [TTj): 

1 £ Q n +\{x) u(z) 
27ri 



(u - I„u)(x) 



i+x(z) z-x 



dz, Vie [-1,1], 



(2-19) 



where Q n+ i(x) = C x +1 (x) or (1 — x )C n _ 1 (x). Consequently, we have 

(u-I n u)'( Xj ) = ±<f Q } +li ^ U{Z) dz, 0<j<n. 

2m h p Qn+l(Z) Z - Xj 



(2.20) 



A crucial component of the error analysis is to obtain a sharp asymptotic estimate of Q n +i{z) 
on E p with large n. This will be the main concern of the forthcoming section. 

3. Asymptotic estimate of Gegenbauer polynomials on £ p 
Much of our analysis relies on the following representation of the Gegenbauer polynomial. 
Lemma 3.1. Let z = \{w + ur 1 ). We have 



g x k g x n _ k w n ~ 2 \ n>0, A > -1/2, 



where 



fc=o 



9a =1, 9 k = 



k + X-1 
k 



r(fc + A) 
*!r(A) : 



Kk<n. 



(3.1) 



(3.2) 



This formula is derived from the three-term recurrence formula (2.1 ) and the mathematical 
induction. Its proof is provided in Appendix [5] 

Remark 3.1. Some consequences of Lemma |3.1| are in order. 



(a) Comparing the coefficients w n on both sides of (3.1|, we find that the leading co- 
efficient of C x is 2 n g x . This can be also verified from (2.1) by the mathematical 
induction. 

(b) If A > 0, then g x > for all < k < n. On the other hand, if A < 0, we find from 
( |2.10| ) and (^ that 



a _ sin(TrA) T(k + A)r(l - A) 



1< k < 



(3.3) 



(c) If A = 1/2, it follows from ((2^ and (|3_2|) that 

(2k)l 



9k 



0<k< 



(fc!) 2 2 2fc ' 

Such a representation for the Legendre polynomial can be found in, e.g., [TT] and 
but the derivation is quite different. 



(3.4) 



(d) If A = 0, then by (2.2 1, 



T n (z) 



2 V 



n > 1. 



(3.5) 
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(e) If A = 1, then g£ = 1 for < k < n. Therefore, the Chebyshev polynomial of the 
second kind has the representation 



U n (z) 



w n+1 — w _ ( n+1 ) 



w — w 

which can be found in |33j . 



w n ^~ 2k = C 1 n {z), n>0, 



(3.6) 



fc=0 



□ 



It is interesting to observe from (3.6| that for A = 1, C^(z)/w n converges to (1 — w 
uniformly for all \w\ > 1, that is, 



E 

fe=0 



w- 2k = 



1 



1 — w 



35, \w\>l. 



In what follows, we show a similar property holds for general A > —1/2 and A ^ 0. More 
precisely, we estimate the upper bound of remainder: 



(i-w- 2 y x - 

9n wn ' 



J2 \9k\P~ 2k ■■= Rn(p,\), (3.7) 



k = l 



k=n+l 



where z G £ p with \w\ = p > 1, and 



< fc = l-%, l<fc< 

9n 



(3.8) 



To obtain a sharp estimate of R n (p,X), it is necessary to understand the behavior of the 
coefficients and {d A fc }jJ =1 , which are summarized in the following two lemmas. 



Lemma 3.2. For A > -1/2, k > 1 and fc + A > 1, 

H "(fc + A) A - 

where C\ = e~™ and C2 = e 151*+*) . 
Proof. Applying the Stirling's formula ( |2.11[ ) to 

(fc + A)r(A).g, A 



T <«(. + ±) 



A\ fc +V2 



(3.9) 



A _ r(fc + A + i) 



r(fc + i) 



leads to (3.9). 



□ 



Lemma 3.3. Let {g£}% =1 and {d* fe }£ =1 be the sequences as defined in ( |3.2[ ) and ( |3.8[ ) 7 

respectively. 

(i) If X> 1, then there holds 

< < i < d, A i!2 < ■ ■ • < < 1. 

(ii) // -1/2 < A < 1 and A ^ 0, iften 

•■■ < < l.9fcl < ••• < l.9il <.9o =1. 

and we have 



< -4.1 < -<2 < ' ' ' < ~<„- 



and 



4,fcfffc I < 1 for 1 < fc < n - 1, n > 3. 



(3.10) 

(3.11) 
(3.12) 
(3.13) 
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Proof. By (3.2), 



0AH 



9 k 



k + \ 

k + l 



for A ^ 0. 



(3.14) 



Thus for A > 1, is strictly increasing with respect to k, which, together with the fact 

g£ > 0, implies 



< d x nk = 1 - < 1, 

9^ 



and 



a n,k+l a n,k ~ 



9n-k 9n-k-\ 
9n 



> 0. 



This completes the proof of (i) . 

The property (3.11) is a direct consequence of (3.14), and (3.12) can be proved in a fashion 
similar to (3.10). It remains to verify (3.13). If k = 1, a direct calculation by using (3.2) 
yields 



<i3i I = |A| 



1 - A 



n — (1 — A) 



< 1, Vn > 3. 



For 2 < k < n - 1, it follows from (|3.2|) and (|3.12| that 

A 

hi 

f) 

This ends the proof. 



fc-2 i _ 1-A 
k-j 



ijA A I (9'n-k i\i A I . 9n—k I A I / TT 1 fc-j 1 M ^ i 

K^fel = (— Y--lJ|5fel<— Y-|Pfel= Ml !_ i4 ! i-A <X 



3=0 n-j 



n-k+1 



□ 



With the above preparation, we are ready to present the main result on the upper bound 
of R n (p,X) in 



Theorem 3.1. Let R n {p, A) with p > 1 be the remainder as defined in (3.7). 
(i) If X > 1, then for all n > m > 1 and 



m+2 ^ (A - 1) ( 2 ib- 1 )' 



have 



Rn{ P , a) < < m (d - P -r x - 1) + ^ {2 \Z/ m P ZT - 

where [A] is the largest integer < A, and 



= _j_ ex ( L_ 

_ r(A) 6XP ll2(m+ 1 



A 



"(A) *A12(m + l + A) 2(m+l), 
(ii) // —1/2 < A < 1 and A ^ 0, then for all n> m> 1 and n > 3, 



i?„(p,A)< |< m | 



2\-A 



-2m 

4— +2p- 2 ». 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



.ffere, i/ie factor d„ m is gwen by (3.8). 
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Proof, (i) For A > 1, we obtain from ( 3.10| ) that 



E d lk9kP~ 2k 



E 9 x kP - 2k 



k=l 
\3~10 



k—m+1 



k=n+l 



E 9 x k p- 2k 



E 9*p~ 2k 



(3.19) 



^,m((l-P 

By Lemma [3^2} 

oo 

E ^P" 2fe < (r(A))- 1 



fc — 1 k—m-\-l 



i)+ E flip 

k—m+1 



fe=n+l 
-2* 



_ A ! ^ (fc + A) 



,,2 A- 



oo 



n 2fc 



k—m+1 

where we used the inequality 1 + x < e x for x > 0. Hence, 

EA -2ft ^ ,i + 1 

fc— m+1 



k— m+1 



/0 



,2 A 



(3.20) 



where A is given by (3.17). One verifies that under the condition ( 3.15[ ), (fc + A) A 1 / p 2k is 
decreasing with respect to fc. Therefore, by (2.121 and (2.13), 



E 

k—m-\-l 



(fc + A) 



A— l r°° 



,2k 



< I (x + Xj^p-^dx 



„2A /-oo 



(21np) ; 



^- 1 e- x dx 



2(m+A) lnp 



^ )X r(A, 2(m + A) In „) < ^^([A] + 1, 2(m + A) lnp) 

[A]!p~ 2m (m + A)*(21np) fc [A]! (m + A)M y> (2 lnp) fc 
\A n - (21np) A I 2 ™ ^ 



[A]! (m + A)M 
(21np) A p2(m-l) • 



tr" 1 ' fc! 

F A:=0 



A combination of the above estimates leads to (3.16). 



(ii) Now, we turn to the proof of the second case: — 1/2 < A < 1 and A / 0. By Lemma 



EH 



R n (p,\) = Y,\dik\\g X k\p- 2k + E Kk\\9 x k \p- 2k + E \9k\p 



2k 



k=l 
i3~ 



k=m+l 



KmlElflfclP" 



2A- 



A;=n+1 



A.-1 



E P~ 2k + KnflnlP" 2 " 
fc— m+1 



-2ft 



A;=n+1 



-2m 

< l< ro ||(l - P^)" A " 1| + ^-r + 2p- 2 ", 
i 1 — 1 

where in the last step, we used the following facts: 

m oo 

E \9 X k \p- 2k < sign(A) Y,9 X kP- 2k = sign(A)((l - p- 2 )- A - l), p > 1, 



k=l 



k=l 
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(note: sign(A) is the sign of A), and \d* „<?^| = \g„ — 1 1 < 2, thanks to (3.11). 



□ 



The estimate in Theorem 3.1 is quite tight and is valid even for small n. By choosing a 
suitable m to balance the two error terms in the upper bound, we are able to derive the 
anticipated asymptotic estimate. 

Theorem 3.2. For any z £ £ p with \w\ = p > 1, and any A > —1/2 and A ^ 0, there exists 
< e < 1/2 such that 



_ 2 \-A C n ( z ) 



< Aip.X)^- 1 +0(n- 1 ), 



where 



^(p,A) = |l-A||(l-p- 2 )- 



1 



(3.21) 
(3.22) 



Proof. We first estimate |<i„ m | in Theorem 3.1 when n — to is large. Using the Stirling's 
formula A2.UI and dO) that 



1 X 

-jn — m 
9n 



= (1 
1 



1 - A 



yH / _ l-Xy-^H , _ to y-ir + (_±\ | 
/ V n — ml V n + A — 1/ L \n — ml> 



V n + A - 1/ I Vn - to/ J 



= 1 



n + 


A - 1 




TO 


n + 


A-l 


(1- 


- A)to 


n 4 


A - 1 


? TO 


- \n e ) 



( 1+ n + A- l + °{n 2 )}{ 1 + °(n-m)}' 



9n- 



9n 



l + (l-\)n e - 1 + o(—^-) 

* < m = (A-i) n £ - 1 + o( n - 1 ). 



One verifies readily that for A > 1 and any < s < 1/2, 



»W 1 

< - 



p2m n 



In n 2 In p 

< f-t 

n £ ~ 1 + e[A] 



(3.23) 



(3.24) 



which, together with (3.161 and (3.231, implies (3.211 with A > 1. 

If -1/2 < A < 1 and A ^ 0, it follows from (feglh that p" 2m < n" 1 for any < e < 1/2. 



This validates the desired estimate. 



□ 



A direct consequence of Theorem 3.2 is that 



lim 



lim y ^ 



k=o y " 



k A 
~9k w 



-2k 



= (l-W~ 



for all z G £ p with |iu| = p > 1, and any A > —1/2 and A ^ 0. 



(3.25) 



Remark 3.2. Based on a completely different argument, Elliott [TH] derived an asymptotic 



expansion for large n near z = 1 (but not near z = —1): C„(z) 



B(n,X) 



, where _B is a series 



( z 2_!)A/2 ; 

involving modified Bessel functions, and some other asymptotic expansions for \z\ large and 
n fixed. Although they are valid for general z off the interval [—1,1], our results in Theorems 
3.1 and 3.2 provide tighter and sharper bounds when z is sitting on £ p . □ 
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As the end of this section, we provide some numerical results to illustrate the tightness of 



the upper bound in (3.21). Denote by 



E n (p;X) 



1 



A(p, A) zgs p 



1 



w 



9n wT 



(3.26) 



To approximate the maximum value, we sample a set of points dense on the ellipse £ p based 
on the conformal mapping z — ~(w + w^ 1 ) of the Fourier points on the circle w — pe 10 . 
We plot in Figure 3.1 in (Matlab) log-log scale of E n (p; A), nT 1 and n e_1 (with e = 0.1) for 
several sets of parameters A and p, and for large n. According to Theorem |3.2| E n should 
be bounded by n e ~ l from above, and it is anticipated to be bounded below by n^ 1 , if the 
estimate is tight. Indeed, we observe from Figure [XT] such a behavior when n is large. 




FIGURE 3.1. E n against n e 1 (with e = 0.1) and n 1 for large n. 



4. Error estimates of interpolation and spectral differentiation 

After collecting all the necessary results, we are ready to estimate exponential convergence 
of interpolation and spectral differentiation of analytic functions. 

Hereafter, the notation a n = b n means that a n /b n — > 1 as n — > oo, for any two sequences 
{o n } and {b n } (with b n ^ 0) of complex numbers. 
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4.1. Gegenbauer-Gauss interpolation and differentiation. We start with the analysis 
of interpolation and spectral differentiation on zeros of the Gegenbauer polynomial C x +1 (x). 

Theorem 4.1. Let u be analytic on and within the ellipse £ p with foci ±1 and p > 1 as defined 
in (2.151, and let (I n u)(x) be the interpolant of u(x) at the set of (n + 1) Gegenbauer-Gauss 
points. 

(i) If A > 0, we have 



max 



i(^)Mi< ^ M :^ . 



r(2A)(p- i) 2 (i + P - 2 )- x P n 



(ii) If -1/2 < A < 0, we have 

max (u — I n u) (x) < 

\x\<l 



cD x \r(\)\M p ^p 2 + p- 2 1 

(p- 1) 2 (1 -p- 2 )- x p n 



(4.1) 



(4.2) 



Here, M p — max zS £ |u(z)|, D\ is defined in (2.7 1, and c = 1 is a generic positive constant. 

Proof. By the formula ( |2.19| ) with Q n+1 = C* +1 and ( |2. 17[ )- ( |2. 18"] ) , we have the bound of the 
point- wise error: 



I (it - I n u)(x)\ < 



\Cn+i( x )\ max ^e£ p \u(z)\ 



2?r 



min ze£p \C* +1 (z)\ J £p \z - x\ 



\dz\ 



< 



< 



M P L(S P ) \C x +1 (x)\ 



2tt(5 



\C x n+1 (z)\ 



(4.3) 



p + p- 



-,s P n +1 (z)\' 



xe [-1,1], n>0. 



Therefore, it is essential to obtain the lower bound of |C* +1 (z)|. Recall that for any two 



complex numbers z\ and Z2, we have \ \zi \ — \z2\1 < \z± — £2!- It follows from Theorem 

— 2 I — A \Cn+li z 



3.2 



that 



which implies 



1 



1 — w 



\9n+l\P n+1 



<A{p, A)n e + 0(n ), 



- A(p, A)?/" 1 - Oin- 1 ) < 



\C r x l+1 (z)\ 

\9n+l\p n+1 



Notice that 
Consequently, 



< 1 - w~ 



l-p- 2 <\l 



i-A 



+ A(p, A)n + 0(n~ ) 



\C' +1 (z)\>c 



n \-l p n+l 



|r(A)| 

where we used (3.9), and the constant c = 1. 



w~ 2 \ < 1 + p~ 2 . 



(l-p- 2 )-\ 



if A > 0, 
if A < 0, 



On the other hand, we derive from (2.5), (2.7) and ( |2.11 1 that if A > 0, 



max|C* +1 (z)|=C£+i(l) = 

\x\<l 



,2A-1 



r(2A)' 



(4.4) 

(4.5) 
(4.6) 

(4.7) 



Hence, a combination of (4.3), (4.6) and (4.7) leads to (4.1). Similarly, for —1/2 < A < 0, we 
use |2?7| to derive E2l). □ 
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Remark 4.1. For A > 0, we obtain from (|2.4|, (|2.8|) and (|2.11|) that 



|C£ +1 (a:)| < 



^(2 + ^) ^ 



r(A) 



(1 - .T 2 )- A / 2 , | Z | < 1. 



(4.8) 



Replacing (4.7) by this bound in the above proof, we can derive the point-wise estimate for 
A > : 

(l-x 2 )- x ' 2 



\(u-I n u)(x)\ < D(p,\)- 



X < 1, 



p" 



(4.9) 



where the positive constant D(p, A) can be worked out as well. It appears to be sharper than 
(4.1 1 at the points which are not too close to the endpoints x = ±1. A similar remark also 
applies to the Gegenbauer-Gauss-Lobatto interpolation to be addressed in a minute. □ 

Now, we turn to the estimate of spectral differentiation. 

Theorem 4.2. Let u be analytic on and within the ellipse £ p with foci ±1 and p > 1 as 
defined in ( |2.15 1, and let (I n u)(x) be the interpolant of u(x) at (n + 1) G eg enbauer- Gauss 

(4.10) 



points {xj} 7 f =0 . Then we have 



,A+2 



max 

0<j<n 



(u- I n u)'{xj)\ < A(p,A)^— 



where the constant 



A(p,A) 



2 C r(A + i)M pV y + P - 2 j(i + p- 2 y 



r(2A + 2)(p- l) 2 



(l-p- 2 ) x 



if A > 0, 
if A < 0, 



(4.11) 



and c, M p are the same as in Theorem 4.1 



Proof. In view of (2.19) and (2.20), it is enough to replace x and C^ +1 (x) by Xj and 
j^C^ +1 (x), respectively, in (4.3). Thus, we have 

m p KC^'ixj)] 



(u - I n u)'(Xj)\ < 



2tt 



\dz\ 



< 



p + p^-2 n 



(4.12) 



< j < n. 



By (2.6) and (4.7), 



max 



|(C A +1 )'Wh2|A||C' A + 1 (l)| 



21 A| 



2A+1 



T(2A + 2) 



which, together with (|4.6|) and (|4.12|), leads to the desired estimate. 

■(C A +1 )'(x) , C A +1 (x)x u{z) 



Remark 4.2. Obviously, by (2.191, 



(u - I n u)'(x) 



1 

2?ri /. 

•'Op 



(z-x) 2 JC^ +1 (z) 



-dz. 



(4.13) 
□ 

(4.14) 



If x =/= Xj, we need to estimate the second term in the summation, which can be done in 
the same fashion in the proof of Theorem |4.1| The first term is actually estimated above. 
Consequently, we have 

n X+2 1 

max (it — I n u)'{x)\ < A(/5, A) h v~ max \(u — I n u)(x) , 

|a;|<l 1 P n Op \x\<l 1 
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where S p is given by (2.18). Hence, by Theorem 



4.1 



/n x ^ 2 \ 

max \ (u — I n u)'(x) \ = 0\ ) 

\x\<l |V ' U| V p n ) 



(4.15) 



In fact, the results for higher-order derivatives can be derived recursively, and it is anticipated 
that 



A+2/c , 



max I (u - 7„u) (fc) (a;)| = 0[ - ), k 

\x\<i |V ; W| V p n J 



> 1. 



A similar remark applies to the Gegenbauer-Gauss-Lobatto case below. 



(4.16) 
□ 



4.2. Gegenbauer-Gauss-Lobatto interpolation and differentiation. We are now in a 

position to estimate the Gegenbauer-Gauss-Lobatto interpolation and spectral differentiation. 
In this case, Qn+i{x) = {l-x 2 )C x +\(x) in pl9)-plQ) The main result is stated as follows. 

Theorem 4.3. Let u be analytic on and within the ellipse £ p with foci ±1 and p > 1 as 
defined in (2.15), and let (I n u)(x) be the interpolant of u(x) at the set 0/(71 +1) Gegenbauer- 
Gauss-Lobatto points. 

(a) We have the interpolation error: 

AcM p ^Jp 2 + p- 2 (l + p- 2 ) x+1 T(A + 1) n A+1 



max (u — I n u) (x) | < 

|x|<l 



(1- p-!) 2 (p- p- 1 ) 2 T(2A + 2) p" 



(b) We have the estimate: 
max \(u — I n uY (xj)\ < 

0<j<n 1 

Here, c = 1 and M p = max ze £ p |m(z)|. 
Proof. For any z € £ p , one verifies that 



8cM pV V +p- 2 (l + p- 2 ) x+1 T(A + 2) n x + 3 
(l-p-i)2(p_ p -i)2 r(2A + 4) p n ' 



(4.17) 



(4.18) 



and 



irrin 1(1-^)6^)1 > ~{p - p^f min \C*±[{z) 



1 



A+l, 



4 



(4.19) 



(a) By j2T9| ) with Q n+ i(a;) = (1 - a^C^ (x) and ( |2T7] -( |2T8| ), we have the bound of 
the point-wise error: 



Uu - L n u)(x)\ < 



\(l-x 2 )C x +l(x)\ 
2tt 



max ze£ \u(z)\ 



\(l-z*)C?l±l(z)\r £p \z-x\ 



\dz\ 



m m p ^/p 2 + p - 2 Mp-p-^c^Kx 



(4.20) 



P- 



un ze£p \C x +\(z)\ 



e [-1,1], n > 0. 



Thus, the estimate (4.17) follows from (4.6) and (4.7). 
(b) Similarly, we have 

M P \[(l-x 2 )C^l(x)Y( Xj )\ 



\(u - I n u)'(xj)\ < 



2tt 



mm ze£ 



(l-z 2 )C^l(z) 



\dz\ 



< 



4M p ^/p 2 +p- 2 



(p-p-^mm zee \C^l(z) 



P + p-'- 
A direct calculation leads to 

[(1 - x 2 )C x t\(x)}' = -2xC x t\{x) + (1 - x 2 )[C x t\{x)}' , 



0<j<n. 



(4.21) 



(4.22) 
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which, together with (4.7) and (4.13), gives 



max [(l-x 2 )C*±l(x)}' * 



M<i 



2(A + 1) 
r(2A + 4)' 



,2A+3 



A combination of (4.6), (4.21) and (4.23) yields the desired estimate. 



(4.23) 

□ 

Remark 4.3. Similar to Remarks |4.1| and |4.2[ we can derive a sharper point-wise estimate and 
analyze higher-order derivatives of interpolation errors. □ 

4.3. Analysis of quadrature errors. Recall the interpolatory Gegenbauer-Gauss-type quad- 
rature formula: 

.1 n . i 

u(x)(l - x 2 ) x - 1/2 dx « u{xt)ujj = j (I n u)(x)(l - x 2 ) x -^ 2 dx, (4.24) 

where the quadrature weights {wj}" =0 are expressed by the Lagrange basis polynomials (see, 
e.g., [MIEO]). Observe that 
-i 

(u - I n u){x){\ - x 2 ) x ~ 1/2 dx <h%max\(u-I n u)(x)\, (4.25) 



|x|<l 



where h$ is given by (2.4). With the aid of interpolation error estimates in Theorem 



4.1 



and 



Theorem 4.3 we are able to derive the quadrature errors immediately. 

4.4. Numerical results. In what follows, we provide two numerical examples to demonstrate 



the sharpness of the estimates established in Theorems 4.2 and 4.3 
4.4.1. Example 1. We take 

«(*) = saVr 



(4.26) 



which has two simple poles at ±i. By (2.16), we are free to choose £ p with p in the range 

1< p < 1 + V2 « 2.414, (4.27) 

such that u is analytic on and within £ p . To compare the (discrete) maximum error of spectral 
differentiation with the upper bound, we sample about 2000 values of p equally from (1, 1 + 
y/2) , and find a tighter upper bound (which is usually attained when p is close to 1 + \ 2 ) . In 
Figure 4.1 (a)-(b), we plot the (discrete) maximum errors of spectral differentiation against 
the upper bounds. We visualize the exponential decay of the errors, and the upper bounds 
and the errors decay at almost the same rate. Moreover, it seems that the bounds are slightly 
sharper in the Gegenbauer-Gauss case. 

4.4.2. Example 2. The estimates indicates that the errors essentially depend on the location 
of singularity (although it affects the constant M p ) rather than the behavior of the singularity. 
To show this, we test the function with poles at ±i of order 2 : 

We plot in Figure [4~T| (c)-(d) the errors and upper bounds as in (a)-(b). Indeed, a similar 
convergence behavior is observed. Indeed, Boyd [5] pointed out that the type of singularity 
might change the rate of convergence by a power of n, but not an exponential function of n. 



Concluding remarks 
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1=1/2 



10 

10" : 
10"' 
10"' 

icr 

10" 
10" 













-•-GG SP err 
















— upper bound for GG 
• - GGL SP err 
- - upper bound for GGL - 








































* 
















N 































































5 10 15 20 25 30 35 40 45 
n 

(c) Upper bound vs Max. -error (Example 2) 
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(b) Upper bound vs Max. -error (Example 1) 
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(d) Upper bound vs Max. -error (Example 2) 



FIGURE 4.1. The (discrete) maximum errors of Gegenbauer-Gauss and 
Gegenbauer-Gauss-Lobatto spectral differentiation against the upper bounds in 
Theorems |4~2| and \H\ with A = 1/2, 3/2. Example 1 (a)-(b), and Example 2 (c)-(d). 
In the legend, GG SP err (resp. GGL SP err) represents the Gegenbauer-Gauss 
(resp. Gegenbauer-Gauss-Lobatto) spectral differentiation error. 



In this paper, the exponential convergence of Gegenbauer interpolation, spectral differenti- 
ation and quadrature of functions analytic on and within a sizable ellipse is analyzed. Sharp 
estimates in the maximum norm with explicit dependence on the important parameters are 
obtained. Illustrative numerical results are provided to support the analysis. For clarity of 
presentation, it is assumed that A is fixed in our analysis, but the dependence of the error on 
this parameter can also be tracked if it is necessary. 
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Appendix A. Exponential convergence of Gegenbauer expansions 
To this end, we discuss the exponential convergence of Gegenbauer polynomial (with fixed 



A) expansions of functions satisfying the assumption of analyticity in [53] , that is, (2.14) in 
the form: 



max 



d k u ^ , /,-! 

dx k 



(x) <C(6) ¥ , 6>1, Vfc>0, (A.l) 



where C(6) is a positive constant depending only on S. We write 

oo „i 
u(x) = ^*}C l x (x) with ut = -r X u(x)Ct(x)(l - xy-^dx, 



h 

1=0 1 



where hf is defined in (|2.4[) 



Proposition A.l. If u[x) is an analytic function on [—1, 1] satisfying the assumption (A.l), 
then for fixed A > —1/2, we have 



max |(tt*u -«)(*) I <0{S,X)-^=-, (A.2) 
|a;|<l 1 1 [2o) n 

where (tt^u)(x) — ^"^qU^C^x), and C(<5, A) is a positive constant depending on 6 and A. 

Proof. The estimate follows directly from Theorem 4.3 in |24j . Indeed, by Theorem 4.3 in 
[24], we have 

i , \ ,, ,, i C(6)T(X+ l/2)r(n + 2A + 1) 

max (?r> - u)(x)\ < A w J: '—^ '-, rc > 1, A > -1/2, 

M<i' nv/A(2i)»r(2A)r(n + A) " 

where A is a generic positive constant independent of A, n and u. For fixed A > —1/2, using 



(|2_llj) leads to 

r(n + 2A + 1) _ n + A + 1 r(n + 2A + 2) 
T(n + A + 1) ~~ n + 2\+ 1 T(n + A + 2) 



A X «+A+3/2 _ A 



e 



< A(n + 2A + 1) A f 1 + -) 

~ v 7 V n + A + 1 / 

< A(n + 2X + l) A e 2 <"+' A + 1 ) , 

where in the last step, we used the inequality 1 + x < e x for all real x. Therefore, we obtain 
the estimate 

\f * w M ^ + A(?t + 2A + 1) A 



This implies (A.2 1. □ 



Remark A.l. The error max| x |<i | (7r A u— | is termed as the regularization error in |24U23j . 
This, together with the so-called truncation error, contributes to the error of the Gegenbauer 
reconstruction. As shown in |24j , an exponential convergence could be recovered from Fourier 
partial sum of a nonperiodic analytic function, when A grows linearly with n. However, the 



estimate (A.2) shows that the regularization error alone decays exponentially for fixed A. □ 
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Remark A. 2. It is also interesting to study the exponential convergence of Gegenbauer ex- 
pansions of analytic functions, which are analytic on and within the Bernstein ellipse £ p in 



(2.15). We refer to [THIIMIES] for the results on Chebyshev expansions. However, the analysis 
for the Legendre expansions is much more involved. Davis [TU] stated an estimate due to K. 
Neumann (see Page 312 of [10]), that is, if u is analytic on and inside £ p with p > 1, then the 
Legendre expansion coefficient satisfies 

limsupKI 1 /" = p~ l . (A.3) 

n— j-oo 

A more informative estimate: 

77 

\u n \<C(p)^, Vn>0, (A.4) 

was presented in the very recent paper |42j . Some discussions on the Gegenbauer expansions 
can be found in [23] . We shall report the analysis for the general Jacobi expansions in [43] . □ 



Appendix B. Proof of Lemma I3TT1 
We carry out the proof by induction. 



Apparently, by (2.1), Cq(z) = 1, so (3.1 )-(3.2 1 holds for n = 0. Similarly, we can verify the 
case with n = 1. 

Assume that the formula holds for C^_ 2 (z) and C^_ 1 (z) with n > 2. It follows from the 



three-term recurrence relation (2.1) that 



k=0 
n-2 

-2fe-2 



-(n + 2\-2)J29k9n-2- k i 

k=0 

71—1 71—1 

= (n + A - 1) 9 X k9n-x- k w n - 2k + (n + A - 1) £ 9^-i-k^ 

n-2 

^(n + 2X-2)Y,9bn-2-k^ n - 



2k-2 



k=0 k=0 
n-2 

-2fc-2 



k=0 



(n + A - l)^-iK + «0 + E D ik 9k9n- k ™ n - 2k , 



where 



k=l 



r,A / i \ i\ ( 9n-k-l i 9k-l | / , r>\ o^Sn-fc-l Sfc-l 

£>n,fc = (n+A-1) — ^ + — - (n + 2A-2)^ ^. 



One verifies from (3.2) that 

9l= n + X ~ l 9n-n D x n . k = n. (B.l) 



This completes the induction. □ 
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